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Abstract 

This paper provides a bridge between the classical tiling theory and the complex 
neighborhood self-assembling situations that exist in practice. 

The neighborhood of a position in the plane is the set of coordinates which 
are considered adjacent to it. This includes classical neighborhoods of size 
four, as well as arbitrarily complex neighborhoods. A generalized tile system 
consists of a set of tiles, a neighborhood, and a relation which dictates which 
are the "admissible" neighboring tiles of a given tile. Thus, in correctly formed 
assemblies, tiles are assigned positions of the plane in accordance to this relation. 

We prove that any validly tiled path defined in a given but arbitrary neigh- 
borhood (a zipper) can be simulated by a simple "ribbon" of microtiles. A 
ribbon is a special kind of polyomino, consisting of a non-self-crossing sequence 
of tiles on the plane, in which successive tiles stick along their adjacent edge. 

Finally, we extend this construction to the case of traditional tilings, proving 
that we can simulate arbitrary-neighborhood tilings by simple-neighborhood 
tilings, while preserving some of their essential properties. 
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1. Introduction 

Because of the constant miniaturization of components, microscopic elements 
in the fields of electronics, engineering or even medicine are becoming more 
and more difficult to construct and to manipulate. A recent approach to work 
around this problem is self-assembly. It consists in "programming" the nano- 
components, so that when starting from an initial pool of selected components, 
they automatically aggregate to form bigger and bigger structures, until they 
eventually form a final complex structure. 
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The first formal models for self-assembly were introduced a decade ago [28[ 
30l l[. I n this framework, self-assembling components were modeled by Wang 



tiles [26|, i.e., unit squares that cannot be rotated and that have "glues" on 
each of their four edges. Two tiles stick to each other if they have the same glue 
on their common edge. By carefully designing the glues, and starting with an 
initial tile called "seed" , complex structures can self-assemble. 

The use of this simple model as a formalization of the process of self-assembly 
allowed the application for theoretical studies of dynamical self-assembly of 
many well-known existing results and techniques concerning "static" tilings [26| 
and cellular automata [16| , such as the undecidable problem of the tiling of the 
plane, and the simulation of a Turing machine by a tile system [il. r20j]. 

Most of the theoretical results on self-assembly presume that each tile in- 
teracts via glues only with tiles in its so-called von Neumann neighborhood, 
which includes the four tiles situated at the North, South, East, and West of 
the tile itself 230 Only relatively few recent results consider more general 
cases, such as larger neighborhood [J], or a three-dimensional neighborhood 
Even the most well-known experimental incarnation of square tiles, the DNA 
tiles [aO, HH, [H, |22L I l-'l , deal only with the von Neumann-sized neighborhood, 



where the DNA single strands located at the corner of each rectangular DNA 
tile allow its interaction with four neighbors only. Other experimental situations 
that could be modeled by self-assembly of tiles, such as atomic or molecular in- 
teractions, potentially include more complex scenarios where the neighborhood 
of a tile is both larger and more complex than the von Neumann neighborhood. 
At the limit, one can consider the case of an arbitrary neighborhood where tiles 
that are not physically adjacent to the main tile may be its neighbors, while 
some adjacent tiles may not. 

In [J], it was proved that, for any directed tile system, any von Neumann- 
neighborhood "zipper" (a tiled rectilinear path) can be simulated by a "ribbon" 
constructed with tiles from a new tile system. A ribbon is a non-self-crossing rec- 
tilinear succession of tiles in the plane, where each tile is required to have glues 
matching with two tiles only: its predecessor and its successor on the ribbon- 
path. The construction that simulated a directed von Neumann-neighborhood 
tiled path by a ribbon, replaced each of the existing tiles by so-called "mo- 
tifs" which traced the contours of the initial tile but where, in addition, bumps 
and matching dents along the motif edges simulated both the matching of the 
glues and the directionality of the path. In other words, geometry of the motifs 
was used to simulate glue matc hing . Note that motifs, as well as ribbons, are 
particular cases of polyominoes [2J,[fij]> i- e -> finite and connected sets of DNA 
tiles. 



1 Total tilings can also be seen as bidimensional subshifts of finite type, and studied as such 
in the field of symbolic dynamics (see [lH for a complete introduction, in dimension 1). In 
this settings, the "neighborhood" is arbitrary, by definition of the subshift. Although some 
results exist for subshifts in dimensions greater than 1 25, 27, 14], they do not apply to partial 
tilings. 
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This polyomino construction led to a conjecture by Jarkko Kari, claiming 
that it is possible to "simulate" an arbitrary-neighborhood zipper by a simple 
two-tile- neighborhood ribbon. A hrst step in this direction was [111 . 1 1 01 ] , wherein 
it was proved that a complex-neighborhood zipper, defined for example on the 
Moore neighborhood (von Neumann plus four diagonal neighbors), can be sim- 
ulated by a ribbon of irregularly shaped tiles, where the shape was used to 
simulate the neighborhood relationship. 

The aim of this paper is to answer the above conjecture positively for the 
case of arbitrary-neighborhood zippers, thus providing a bridge between the 
classical work in tiling theory or cellular automata and the realistic complex- 
neighborhood self-assemblies that exist in practice. We namely prove in Corol- 
lary [3T0] that for any neighborhood, zippers can be simulated by simple poly- 
ominoes connected to each other end-to-end to form a ribbon that essentially 
traces the same path. The main idea used in our simulation is that each exist- 
ing tile can be replaced by a polyomino, where the shape of the polyomino is 
used to simulate the communication between a tile and its adjacent or distant 
neighbors. We also show that, by the design of the shapes of the polyominoes, 
we can transmit information at a distance, sometimes across other information 
pathways, without violating the non-self-crossing feature of the ribbon. Such 
situations where information pathways cross are inherent in, for example, Moore 
neighborhoods where, e.g., the communication channel between a tile and its 
Northeast neighbor "crosses" the communication channel between its North 
neighbor and its East neighbor. 

We also explain how the simulation of zippers by ribbons can be modi- 
fied to simulate arbitrary-neighborhood tilings by von Neumann-neighborhood 
tilings. The idea is to modify the polyominoes, adding new constraints so that 
the two-tile neighborhood can be replaced by a von Neumann neighborhood 
(Corollary 14.21) . The main significance of this simulation, as opposed to more 
intuitive ones where some "supertiles" are used to transfer information, is that 
it applies to all tilings, even when they are partial. Besides, some essential 
properties of the initial tiling are preserved by the simulation. We prove that 
this is the case for partial and periodic tilings, and in a more restricted setting, 
for convex tilings. 

The paper is organized as follows. In Sect. [21 we recall basic definitions and 
give a formal definition of a simulation. Then, Sect. [3] describes our construction 
in the case of zippers, starting with the general idea of simulating a zipper by 
a ribbon, by sketching the proof from Q , in the simple case of a von Neumann 
neighborhood. We also highlight the technical difficulties related to crossing 
of information pathways, that prevent this technique from being transferable 



without modifications to the case of the Moore neighborhood Then, we 

prove our result for the case of arbitrary linear neighborhoods. This construction 
is eventually generalized to arbitrary neighborhoods to prove the main result 
of this section. Finally, in Sect. 01 we adapt the construction to the simulation 
of regular tilings, and we study how partial, periodic, line or column convex 
properties are preserved. 
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2. Definitions 



First, we give some basic definitions on tilings, and then we introduce more 
technical definitions to describe the principles of simulation by polyominoes. 
We finally illustrate this on a concrete example, the simulation of total tilings. 

2.1. Tilings, Ribbons, and Zippers 

Historically, Wang tiles [26[ were defined as oriented unit squares, on the 
border of which were 4 glues (colors) used to stick them to neighboring tiles, 
provided the glues matched. We generalize this notion to arbitrary neighbor- 
hoods [3,13]; i-e., neighborhoods which can contain other tiles than the North, 
South, West and East tiles. In this paper, a neighborhood TV C Z 2 is the set of 
relative coordinates of the neighboring tiles, such that 

• N is finite; 

• (0, 0) g" N (a tile can not be one of its neighbors); 

• G N => (— i, —j) 6 N (if a tile t' is a neighbor of a tile t, then t has 
to be a neighbor of t'). 

For example, the usual 4-tile neighborhood (called von Neumann neighbor- 
hood [3) is N vN = {(0,1), (1,0), (0,-1), (-1,0)}; the 8-tile Moore neighbor- 
hood is N M = {(0, 1), (1, 1), (1,0), (1, -1), (0, -1), (-1, -1), (-1, 0), (-1, 1)}. 

For the following definitions, a neighborhood N and a finite set of glues X 
are fixed. A tile is a tuple t — j)gjv wnere each ti_j G X. Intuitively, 

tij is the glue that will be used to match the neighbor located at position 
relative to the tile t. A tile system T is a finite set of tiles, used to build 
larger structures. A tile t £ T sticks at position to a tile t' G T if the 

corresponding glues match, i.e., Uj — 

A (T, N)-tiling using tile system T in neighborhood N (or N -neighborhood 
tiling, or simply tiling when the tile system and its neighborhood are known 
without ambiguity) is a mapping r : D — > T, where D C Z 2 is a subset of 
the plane, which associates every position (x, y) G D with a tile, such that all 
tiles stick to their neighbors, i.e., for all (x,y) € D, for all G N such that 
(x + i,y + j) G D, t(x, y)ij = t(x + i,y + Note that tilings are often 

referred to as "valid" tilings in the literature. The set of all (T, 7V)-tilings is 
denoted It.n- 

When D = Z 2 , the tiling is said to be total, otherwise it is partial. In 
addition, if D is finite then the partial tiling is also finite. A tiling r : D — > T is 



connected if its domain D is 4-connected. As suggested in 24|, we call polyomino 
a finite and connected tiling. 

According to the usual definition, a total tiling r is periodic if it admits 
a horizontal and a vertical period Ph,Pv € N+, i.e., for all x,y G Z, r(x,y) = 
T(x+ph, y) = t{x, y+p v ). We extend this definition to partial tilings as follows: 



2 The usual definition of a polyomino (see for example Q) refers to a domain D C Z 2 , 
while here we add to this notion a mapping to tiles. Besides, polyominoes are often defined 
as simply connected sets of squares, i.e., without holes, which is not required in this paper. 
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a tiling r : D — > T is periodic if it admits a horizontal and a vertical period 
Ph,Pv £ Z such that for all (x, y) £ D and for all a,/3 € Z, (x + ap/j, V+fiPv) € P 
implies t(x, y) = T{x+aph,y+j3p v ). Note that with this extended definition, all 
finite tilings are periodic (the periods should be "larger" than the tiling itself). 

A tiling t : D — > T is line convex [resp., column convex] if (x±,y), (x2,y) G D 
and x\ < X2 [resp., (x,yi), (x, 2/2) S -D and yi < j/2] imply that for all x\ < x < 
X2 [resp., for all y\ < y < 2/2], (x, y) £ D. A tiling is convex if it is both line 
and column convex. 

Two positions of the plane u, v € Z 2 (or, by extension, tiles of a tiling) are 
said to be adjacent when they are neighbors in the von Neumann sense, i.e., 
v — u £ N v n ■ A path is a sequence of adjacent positions of the plane. Formally, 
a path is a function P : I — > Z 2 , where / C Z is a set of consecutive integers, 
such that for all 1 e /, P(i) and P(i + 1) arc adjacent. For a given tile 

system T, a T -tiled path using tile system T is a sequence of adjacent tiles from 
T, i.e., a pair (P,r) where P is a path and r : range(P) — > T a mapping from 
positions to tiles. We say that a T-tiled path is finite when dom(P) is finite, 
and we may also call a finite T-tiled path a polyomino since it is connected (by 
definition of the path). 

For a tile system T defined in von Neumann neighborhood N v n, a T-tiled 
path (P, r) is a T -ribbon using tile system T (simply called ribbon when T 
is known without ambiguity) if P is injective (non-self-crossing) and for all 
dom(P), r(P(i)) v = r(P(i + 1))_„, with u = P(i + 1) - P(i) e jV vJ v. 
The glue r(P{i)) v is called the output glue of r(P(j)), while r(P(i + 1))_« is the 
input glue of r(P(i + 1)). Informally, a ribbon is a sequence of adjacent tiles 
which stick to their predecessor and successor only (see Fig.Q ^a)] ). Consequently, 
r is not necessarily a tiling. 

A T-tiled path (P, r) is a (T, N)- zipper using tile system T in neighborhood 
N (or N -neighborhood zipper, or zipper when T and N are known) if P is 
injective and r is a (T, A^)-tiling. A zipper can be seen as a tiling with an 
additional notion of unique input and output for every tile (except the first 
which has only an output and the last which has only an input), each input 
being connected to the output of an adjacent tile (see Fig. [ ^b)[ ). Note that 
zippers can be defined in arbitrary neighborhoods, since it is not required that 
adjacent glues match. The set of T-ribbons is denoted £Ht, the set of (T,N)- 
zippers 3t,at- 

2.2. Poly-Tilings and Simulations 

In this section, we give the main idea of the simulation of a tiling by another, 
as well as for the simulation of a zipper by a ribbon. Informally, a simulation 
is a function which associates each tiling [resp., zipper] using the "initial" tile 
system in the "initial" neighborhood with a tiling [resp., ribbon] using a "new" 
tile system in a "new" neighborhood. The simulation must also allow to recover 
the initial object without ambiguity. To build this new object, we first associate 
every initial tile with unique polyominoes. Then, we catenate these polyominoes 
to form the new, bigger, object, called poly-tiling or poly-ribbon. 
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(a) A ribbon. (b) A zipper. 



Figure 1: A ribbon, and a von Neumann-neighborhood zipper. The underlying path is drawn 
in gray. The only difference occurs at the bottom-right where glues / and n are not required 
to match in the ribbon, while in the zipper they do. 



In the sequel, for a set D 6 Z 2 and a vector u, we denote D + u = 
{(x, y) + u | (x, y) 6 D} the set D translated by u. 

Definition 2.1. A poly-tiling t of scale s, using tile system T in neighborhood 
N, is a mapping r : D —> 1t,n, with BcZ 2 , such that 

(i) for all r(u) is a finite (T, iV)-tiling; 

(ii) for all u,v E D, u ^ v, [dom(Y(u)) + szt] (~l [dom(r(w)) + si;] = 0; 

(Hi) for all u,u' £ D, v E dom(r(it)) and v' € dom(r(«')), let w = (su' + v') — 
(su + v); then w G N implies t(u)(v) w = r(u')(w')_ u ,. 

Let us discuss the items in this definition from an informal point of view, 
as shown in Fig. [5] We refer only to the particular case in which all t(u) are 
connected and thus polyominoes, since this is what we will construct. In this 
case, the definition implies that a poly-tiling of scale s is 



(?) a juxtaposition of polyominoes on a grid of size s; 

(ii) such that all these polyominoes do not overlap once shifted to their actual 
position on the grid, which is achieved by shifting t(x, y) by sx units to 
the right and sy units upwards; 



Hi) neighboring polyominoes stick, i.e., if two tiles of two different polyominoes 
become neighbors after being shifted, they have to stick. 
These characteristics allow to consider them as usual "valid" tilings without 
overlaps, as explained in the following remark. 

Remark 2.1. A poly-tiling r of scale s can be easily transformed into a "regular" 
tiling t', as depicted in Fig. [5J according to the following formula: 

r'(i,j) = (T(x, y ))(i',j') , 

for any x,y,i,j,i',j' € Z which verify i = sx + i' ', j = sy + j' and (i',f) G 



dom(r(a;, y)). Because of condition (ii) in the definition of poly-tilings, given 



any pair (i,j), if x,y,i',j exist then they are unique. Because of (i) and (Hi) 
t' is indeed a tiling since all glues match. 
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Figure 2: Poly-tiling t of scale s = 4 seen as a tiling. The dashed line surrounds the tiling 
t(x, y). Note that some tiles from r(x, y) may be located outside of the "box" of size s situated 
at (sx,sy), and conversely that some tiles of this box (here, t 12 ) may not belong to r(x,y). 



We can slightly modify this notion to define poly-ribbons. Basically, we 



replace the property of sticking at the macro-level (Hi) by a notion of macro- 
path. 

Definition 2.2. A poly-ribbon of scale s, using tile system T (in neighborhood 
N v n), is a pair (P,r) where P is an injective path and r : range(P) — > Dty is 
such that 

(i) for all u 6 range(P), r(u) is a finite T-ribbon; 

(ii) for all u, v S range(P), u ^ v, [dom(r(u)) + su] (~l [dom(r(w)) + sw] = 

(where we define dom(r(u)) = dom(r'), if r(u) = (P',r')); 
(Hi) for all i,i + l e dom(P), let r(P(i)) = (Pi,ri) and r(P(i + 1)) = (P 2 ,r 2 ) 
be two consecutive ribbons. Let u\ = Pi(max(dom(Pi))) be the last 
position of Pi, 1*2 = P2(min(dom(P 2 ))) be the first position of P 2 , and 
v = [U2 + sP(i + 1)] — [u% + sP(i)] be their relative position. Then, we 
must have v € N v n and ri(ui) v — r 2 (u 2 )-t,. 



The new requirement (Hi) enforces that the ribbon follows a path at the 
macro-level, "jumping" from one ribbon to the other. Therefore, in a similar 
way as what was done in Remark l2.ll we can see a poly-ribbon as a ribbon. This 
explains why in our constructions we build poly-tilings [resp., poly-ribbons] for 
simplicity, they can in turn be considered as particular tilings [resp., ribbons] 
from the set 1t,n [resp., SHt]- 

Let us now formalize the intuition that a simulation is a function ip which as- 
sociates tilings with unique poly-tilings [resp., zippers with unique poly-ribbons], 
by means of a function if transforming every tile into a set of unique polyomi- 
noes. These polyominoes will be aligned on a grid of size s, each of them re- 
placing a tile from the initial object at the same position, to form a poly-tiling 
[resp., poly-ribbon] corresponding to the initial tiling [resp., zipper]. 

In the following, for a set S, we define V(S) = 2 s as the set of subsets of S. 
Let tp : T — >• V(1t',n<) be a function which associates a tile t with a set of finite 
(T', 7V')-tilings of scale s, such that for all t G T, for all r G <p(t), (0, 0) € dom(r) 
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(the tile r(0, 0) is called the reference tile), and such that if t ^ t' , for all 
r G ip(t), t' e f{t'), u G dom(r), and it' G dom(r'), then r(u) 7^ t'(u') (all 
tiles from all tilings in ip(t) are specific to tile t). In our constructions, the 
finite tilings tp{x, y) will be connected polyominoes. Let us introduce one more 
mapping Xv '■ ^T',N' — > %T',N', which associates any (T', 7V')-tiling r with a 
poly-tiling x v (t), by removing incomplete or misaligned polyominoes as follows. 

1. First, discard all the tiles from t which do not belong to a complete 
polyomino <p>(t), for some t € T. We obtain n which is a restriction of r. 

2. Then, discard all the polyominoes in t\ which do not have their reference 
tile positioned at (sx, sy) for some 1,1/gZ (this is to avoid misalignment 
in non-connected tilings), to obtain Xy>( r )- 

Note that x<pi T ) is the empty tiling for infinitely many tilings (for example, all 
the ones for which the reference tiles are badly positioned). 

Definition 2.3. Given such functions (p and Xtp> & tiling-simulation of a tile 
system T by a tile system T' is a mapping ip : %t,n — > 'P(x<p(^T',N'))> such 
that ip associates any (T, iV)-tiling r : D — » T with a non-empty subset of poly- 
tilings (hence tilings) from {r' : D — > %t',N' I T '{ x ,y) S <p( T ( x > with the 
additional constraint that for all (T', iV')-tilings r', there exists some unique 
valid (T, A/>tUing r such that X v i T ') e ^( r )- 

Note that the constraint on ip saying that tiles are specific implies the in- 
jectivity of ip and of ip. In fact, we even have the following stronger result: 
ip(r) n ip(T') = for two different (T, iV)-tilings r and t' . 

Let us now see how to recover an initial tiling from any (T", iV')-tiling, once 
a tiling-simulation ip is defined. Given a poly-tiling r' € iP{t), one can restore 
unambiguously the initial (valid) tiling r. Indeed, if T f (x,y) denotes the poly- 
omino whose reference tile is located at (sx,sy), then r(x,y) = </? _1 (S'), where 
S = ip(r(x, y)) is the only image set of ip containing t'(x, y). This tiling r must 
be valid (all glues matching) in order for ip to be a simulation, which is enforced 
by the last constraint of Definition 12.31 (since in this case, Xip{ T ') = r ')- 

Moreover, for any other tiling r' € %t',N' such that there is no tiling r 
with r' G iP(t), the function x^ allows to recover from this situation: as in the 
previous case, by construction of (p, there exists a unique (T, 7V)-tiling r such 
that X<p{ T ') £ ip(r). This tiling r can be recovered as explained in the previous 
paragraph. From the last constraint in Definition 12.31 we also deduce that r is 
valid. 

The simulation of a zipper by a ribbon is defined similarly. Let (p : T — > 
V(9Kt') associate a tile t with a set of unique finite T'-ribbons, , such that for 
all t e T, for all (P, r) G (p(t), (0,0) £ dom(r), and such that if t ^ t' , for all 
(P,r) G tp(t), [P',r') £ <p(t'), u G dom(r), and u' G dom(r'), then r(u) ^ r'(u'). 
Also let Xh> '■ 9^T' — > 5Ht' , which associates any T'-ribbon with a poly-ribbon as 
previously, by removing incomplete and misaligned polyominoes. 

Definition 2.4. A zipper- simulation of the tile system T by the tile system 
T" is a mapping ip : 3t,n — > ~P{Xv(^T')), such that ip associates any (T,N)- 
zipper (P, r) where r : range(P) — > T with a non-empty subset of poly-ribbons 
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(hence ribbons) from {(P, r') \ r' : rangc(P) — > D\t', r'(x, y) £ ip(r(x,y))}; with 
the additional constraint that for all T'-ribbons (P', r'), there exists some unique 
valid (T, iV)-zipper (P,r) such that X<p{( P '' r ')) e i>(( p >r))- 

The recovery procedure proceeds as previously, by turning a T'-ribbon into 
a poly-ribbon with x^, and recovering the zipper tiles with <p~ l . 

Provided a tiling-simulation -0 exists, we say that t' simulates t if Xv>[ r ') € 
t/>(t). Similarly, for a zipper-simulation ip, the ribbon (P',r') simulates the 
zipper (P, r) if Xp(( p '> r')) G ^((P, r)). Finally, when no ambiguity is possible, 
we simply use the term simulation to refer to the simulation of tilings or zippers. 

2.3. Total Tilings Constructions 

In 0] , the authors use a construction which can be adapted to the simulation 
of total tilings. Their idea is to replace each arbitrary-neighborhood tiles by von 
Neumann-neighborhood tiles, these new tiles (called supertiles to avoid misun- 
derstandings) being groupings of several of the initial tiles, such that supertiles 
group all the tiles which belong to the neighborhood of the initial tile. This 
leads to the following statement. 

Proposition 2.1. Let T be a tile system in arbitrary neighborhood N , with set 
of glues X . There exist a simulation ip and a tile system T' in neighborhood 
N v n, with set of glues X' , such that any total (T, N)-tiling can be simulated by 
a total (T" , N v N)-tiling. 

Proof. Let r be a (T, iV)-tiling. We replace tiles from T by new tiles from 
T' = T^ N \ +1 called supertiles, each of them encoding all the tiles from TV and 
the tile itself. The 4 von Neumann glues encode the neighbors located at N n 
i > 0} (East), Nn{(i,j) \ i<0} (West), Nn{(i,j) | j > 0} (North), and 
N(l{(i,j) I j < 0} (South). Therefore, two supertiles stick if the matching glues 
encode the same tiles, allowing to transmit information to distant neighbors. 
Figure [3] illustrates this construction in the case of Moore neighborhood Nm- 
In this figure, we use the following notation: given a position (x, y) G Z 2 , let tile 
t = t{x iV ) e T. The symbol N(t) [resp., S(t), E(t), W{t)] represents the tile 
t(x, y + 1) [resp., t(x, y — 1), t{x + 1, y), t(x — 1, y)]; and denote XY = X o Y 
the composition of any of the functions X,Y £ {N 7 5, E, W}. 

In this case, the polyominoes simulating the original tiles are made of only 
one supertile. Besides, for each tile t £ T, the set tp(t) contains all the one-tile 
tilings consisting of a supertile centered on t, with all the admissible neighbors 
of t. The function x<p is the identity, since all polyominoes are made of one 
supertile, they are necessarily complete and aligned. Any total tiling r £ 1t,n 
is then simulated by an appropriate poly-tiling i/>{t), seen as a total (T',7V„jv)- 
tiling t' £ %t',n vN - 

The converse is trivial, any total (T', A^jv)-tiling can be translated back into 
a unique valid total (T, iV)-tiling by selecting only the central symbol. □ 
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Remark that this technique can not be applied directly to the simulation of 
total zippers. Indeed, one would need to be able to transfer an arbitrary amount 
of information when simulating zippers such as the one represented in Fig. 2J 
since this amount would depend on the length of the zipper before it goes back. 



Figure 4: Von Neumann-neighborhood zipper, going backwards. The longer it is, the more 
information (represented by the circles) would need to be encoded in the supertiles. 



Also note that this simulation does not work for partial tilings, since in the 
absence of intermediate supertiles the information can not be transmitted. For 
example, in Fig. [3l if the top-right and bottom-left supertiles are omitted, it is 
not possible to ensure that the bottom-right tile is centered on SE(t). 

This issue can be partially solved [29], by introducing a new symbol "no 
tile" in the supertiles, when a tile is absent. In this case, there will be supertiles 
at each position of the plane (possibly consisting only of symbols "no tile"). 
However, this would cause every partial tiling to be simulated by a total tiling, 
which is not suitable in most applications. 
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3. Simulating Arbitrary-Neighborhood Zippers by Ribbons 

Here, we prove that there exists a simulation such that any zipper, using a 
tile system in arbitrary neighborhood, can be simulated by a ribbon using an 
appropriate tile system in von Neumann neighborhood. The objective of our 
constructions will be to dchnc the function ip, in such a way that the polyominocs 
it produces are unique and do not overlap if and only if the initial tiles of the 
zipper stick. 

The final construction being quite complex, we introduce the technical dif- 
ficulties progressively. First, we recall known results which present the basic 
principles of the simulations, and solve the problem of crossings. Then, we deal 
with linear neighborhoods (all neighbors are on the same line), and finally we 
prove the general result in Corollary 13. 101 

3.1. Preliminary Results 

First we recall basic results of simulations of zippers. In 4], the authors 
prove a fundamental result on the simulation of zippers by ribbons. They de- 
scribe a method used to simulate bi-infinite zippers using "directed" tiles in von 
Neumann neighborhood by ribbons. The result can be extended to arbitrarily 
long zippers and standard tiles, as recalled here. 

Theorem 3.1 ((Hi). Let T be a tile system in neighborhood N v n, with set of 
glues X . There exist a simulation tp and a tile system T M in neighborhood N v n, 
with set of glues X^, such that any (T,N v n)- zipper can be simulated by a T M - 
ribbon. 

Proof. The key of the proof is the construction of the simulation <p which as- 
sociates each tile with polyominoes. The basic idea behind these polyominocs 
is to replace every tile from T by a unique shape (Figs. E ^a)] and | jfb")j ). Glues 
are replaced by bumps (the tile is raised) and dents (the tile is dug), a different 
glue leading to a different bump or dent. A way to uniquely code the glues is to 
change the vertical or horizontal position of the bump or of the dent, depending 
on the glue. Then, like in a jigsaw puzzle, two shapes stick if their adjacent 
bump and dent fit into each other. Therefore, a "ribbon" of these shapes would 
simulate a (T, ./V^^-zipper, since the bumps and dents imply that the sides 
unconstrained by the ribbon have to match. 

The second step of the construction of ip is the definition of a new tile system 
Tfj,, in von Neumann neighborhood and with a new set of glues, which will be 
used to build the polyominoes simulating the initial tiles from T. This leads to a 
path P, starting from the middle of the side corresponding to the input direction, 
and leaving at the middle of the side given by the output direction (Fig. O ffc)] ). 
This path draws the contour of the shape, including bumps and dents, while 
the central part of the path is used to reconnect the input and output sides. Its 
position is determined by the point at the Southwest corner for example, which 
we choose as the reference tile located at (0, 0). The path is "filled" by new tiles 
from the set T M (we call these microtiles to avoid misunderstandings) using a 
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mapping r : range(P) — > T^. The finite T^-tiled path (P, r) is the polyomino 
associated with the initial tile. Remark that for one initial tile, there can be 
12 different polyominoes, corresponding to the same shape but with 4 possible 
input and 3 possible output directions. 




(a) Tile of the original zipper (b) Coding of the glues into (c) Polyomino used to simulate 
with its four glues. bumps and dents. the original tile. 

Figure 5: Simulation of a von Neumann-neighborhood zipper using a ribbon of microtiles. 
The three steps of the simulation of a single tile are represented, in each one the gray arrow 
indicates the direction of the underlying path of the zipper. 



The polyomino should be a ribbon, so we should give some details about 
its construction. The first microtile (called the input microtile) has its input 
side colored with glue (g,d), where g is the input glue of the initial tile and d 
the direction of the path (West-to-East, etc.); similarly the output side of the 
output microtile encodes the glue of the output side of the initial tile and the 
direction. For all other microtiles, the input matches the output of the previous 
tile, so that our polyomino (P, r) is a T^-ribbon. The input and output glues 
are unique among all the polyominoes, so this polyomino is the only possible 
ribbon using the tile system T M , once the input microtile is given. Besides, to 
ensure that no interference occurs, the two sides which are not colored yet have 
a glue that matches nothing (for example glue null\ on the West or North sides, 
null?, on the East or South sides). 

For a tile t € T, the set cp(t) contains the 12 polyominoes described above. 
Then, obviously, if a T-tiled path is a (T, iV„jv)-zipper, one can find a "unique" 
poly-ribbon consisting of the catenation of polyominoes constructed by ip, which 
is a T^-ribbon. It is not really unique, since at the extremities of a finite zip- 
per, 3 different polyominoes corresponding to the 3 possible input (or output) 
directions are admissible. 

Conversely, because of the careful design of the glues from X^, any T^-ribbon 
can be seen as a poly-ribbon, since the glues forming the polyominoes appear 
only once and guarantee that only polyominoes can be formed. A X^-ribbon 
may have up to two incomplete polyominoes at the extremities, but if we cut 
them off using x<p > then it can be represented as a unique poly- ribbon. Then, if 
a poly-ribbon using tile system exists, Definition 12.21 implies that 

• all polyominoes do not overlap (condition (it) I, hence that the glues they 
simulate match on the four sides; 



the polyominoes stick on their input/output tiles (condition (Hi)), hence 
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that the polyominoes follow a path. 
Using this path and these glues, one can uniquely restore the initial (T,N v n)- 
zipper. □ 



The following remark states an important property of our construction. In 
fact, the simulation ifj we just constructed can be considered as bijective. 

Remark 3.1. Let 1Z be the equivalence relation which states that two T^-ribbons 
are equivalent if they represent the same (T, iV„jv)-zipper. Then, the simulation 
ip is a bijcction between the set of (T, iV„jv)-zippers ~5t,n vN and the set of T^- 
ribbons quoticntcd by 1Z, SHt^/tc.- 

In addition, the polyominoes used in the above simulation are rectilinear 
polyominoes, i.e., a simple sequence of tiles outlining a shape. These are a 
particular case of general polyominoes, making the simulation more powerful. 

Remark 3.2. Theorem 13.11 allows to transfer information on a ribbon in all 4 
directions. This result, in conjunction with Proposition 12.11 allows a simple 
construction in the case of arbitrary-neighborhood total zippers, by first reduc- 
ing to a von Neumann-neighborhood total zipper with the supertiles technique, 
and then simulating it by a two-tile-neighborhood ribbon. 

The simulation of a zipper in Moore neighborhood Nm (i.e., adding diagonal 
communications) is more complex, because diagonal glues cross, as illustrated 
in Fig. El 



NW(ty 


- N(t) - 




— 1—5 

W(t) - 


( 1 ) 
- t - 


< 1 — 


1 ) 


C 1 > 


< 1 



Figure 6: Communication between tiles in Moore neighborhood. 

The consequence is that diagonal bumps would also have to cross each other. 
This issue is solved in [ll| using a method that we summarize here, because it 
will turn out to be useful in our constructions. 

Theorem 3.2 Let T be a tile system in neighborhood Nm, with set of 

glues X . There exist a simulation ip and a tile system T M in neighborhood N v n , 
with set of glues X^, such that any (T, Nm)- zipper can be simulated by a T M - 
ribbon. 
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Proof. The global idea of the simulation is the same as for the proof of Theo- 
rem l3.ll we need to define a function ip for all tiles from T. In Figs . Qfa)] and E flfbJI 
we present a picture of the shape that can be used to simulate an initial tile of 



T. It differs slightly from the shape described in [111], but the idea is the same 



and this new shape is an introduction to our results. 



(-1,1) (0,1) 



(i.i) 



(-1,0) 



(1,0) 



(-1,-1) (0,-1) (1,1-1) ' — ' ' — ' ' — ' 

(a) A single shape. Next to the communication places (b) The shape and its 8 neighbors, 
is written the neighbor it communicates with. 



Figure 7: Shape used to simulate a tile in Moore neighborhood, filled with gray for clarity 
sake. The communication bumps and dents are represented by short thick lines. 

This shape is turned into a polyomino using new microtiles from T M , taking 
into account the input and output directions. Then, all but one of the commu- 
nications can be done as previously, by modifying the position of bumps and 
dents to simulate different glues. The only issue is the communications between 
the Northwest and Southeast neighbors (circled on Fig. [ ffib)| . We have to be 
able to relay information about these diagonal glues without the physical touch 
between edges of tiles. According to the notation from Fig. HI we need to check 
the glue compatibility between the central tile t and its Northwest neighbor 
NW(t), as well as between the tiles N(t) and W(t). This can be accomplished 
by a geometrical construction such as the one in Figs. | ^a)| and[ ^b)| 



NW(t) 
I I I I 



NW{t) 
I I I I I 



II 



II 



W(t) 



N(t) 



t t 

(a) Glues matching between t and NW(t). (b) Glues not matching: the ribbons overlap. 

Figure 8: "Crossing" of information. The top ribbon segment is part of NW(t), the bottom 
ribbon segment is part of t, and the space in between can be filled by 2 layers of microtiles 
for W(t) to actually reach N(t). 
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This construction checks the match between the W(t) and N(t) in the old- 
fashioned way, by physically matching the adjacent bump and dent of the corre- 
sponding polyominoes. The novelty of this construction is that the glue-match 
between t and NW(t) is accomplished without the respective polyominoes ever 
touching. Moreover, this construction works even in the case of partial tilings 
where the tile Wit) might be missing. This is accomplished by carefully de- 
signing the shape and space between the bump and the dent so that, whether 
or not the tile between these polyominoes is present, their shapes will be com- 
patible and not overlap if and only if the glues of the corresponding tiles were 
compatible. In order to achieve this and cross the 2 layers of microtiles forming 
the polyomino Wit), the bump of t should be 3 microtiles high, and the dent 
of NW(t) should be 8 microtiles wide and 3 microtiles deep. In the sequel, the 
inner layers of Wit) are called bridges. 

For any tile t e T, if a polyomino in the set (pit) has a bridge, then (pit) 
should contain the same polyomino with the bridge at all possible locations, 
matching all possible glues. Indeed, the bridges do not participate in a com- 
munication, they are just a kind of "relay" which should be able to match any 
glue. Therefore, (pit) contains \X\ copies of each of the 12 polyominoes given 
by the above construction, for the \X\ possible locations of the bridges (one per 
glue). 

The end of the proof, which consists in restoring a (T, ./Vm )-zipper from any 
T^-ribbon, is similar to the proof of Theorem 13.11 and is left to the reader. □ 

Remark 3.3. A simpler construction [2j using "pitcher-tiles" (see Fig. ^ solves 
the problem of information crossing when simulating Moore-neighborhood zip- 
pers by ribbons, but only when the zipper is a total tiling. In that case, one 
could use (with the notation from Fig. [9]) the spike of the Wit) polyomino 
which conveys information to N(t) as "information carrier" to transmit infor- 
mation from t to NWit). This construction does not work as such in the case 
of zippers which are partial tilings, because the "carrier" Wit) tile might be 
altogether absent. 




Figure 9: "Pitcher-tiles" used to relay diagonal information, for total-tiling zippers only. 

This idea is similar to the one suggested in Remark 13.21 since it uses in- 
termediate polyominoes to relay information. The main difference is that the 
polyomino-based simulation is more general and can be used for arbitrary neigh- 
borhoods, while these pitcher-tiles were designed specifically for the Moore 
neighborhood. 
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3.2. Simulating Arbitrary Linear Neighborhoods 

We call linear neighborhood a neighborhood N such that if S N then 
j = 0. Although this is a sub-case of the general case studied in the next section, 
we will explain it in detail since it lays the base of the general study, and it is 
much simpler to describe and understand. Also note that this result was already 
announced in [10( , but the polyomino used there was not easily generalizable to 
arbitrary neighborhoods. First, we state an initial remark which allows us to 
consider only connected linear neig hborhoods N n = {(i,0) £ Z 2 | < \i\ < n). 

Remark 3.4- For a given set of glues X, any tile system T defined in linear 
neighborhood N can be replaced by an equivalent tile system T' in an appro- 
priate connected linear neighborhood N n . Indeed, let n be such that N C N n , 
let g G X be an arbitrary "dummy" glue and T' — Cg{T) a tile system in 
neighborhood N n , where £ s : T — > T" is defined for all t S T and S N n by 



Then, clearly, r : D T is a (T, AO-tiling if and only if r' : D ->■ T' defined by 
r'{x,y) = Cg(r(x,y)) is a (T', A„)-tiling. 



First, we state a lemma which generalizes to an arbitrary number of inner 
layers the crossing operation detailed in the proof of Theorem 13.21 The gadget 
introduced in this lemma will be helpful in the next constructions. 

Lemma 3.3. In order to fit a bump and a dent spaced by k layers, the bump 
must be k + 1 microtiles high, the dent must be 2k + 4 microtiles wide and k + 1 
microtiles deep. 

Proof. The result is obtained by an immediate generalization of Fig. [ ^a)| to k 
white layers instead of 2. □ 

We now prove the main result of this section for neighborhoods N n . 

Theorem 3.4. Let T be a tile system in connected linear neighborhood N n , with 
set of glues X . There exist a simulation tp and a tile system T) t in neighborhood 
N v n, with set of glues X^, such that any (T, N n )-zipper can be simulated by a 
T^-ribbon. 

Proof. As previously, we are going to replace a tile of the (T, A ra )-zipper by a 
shape using a function ip, leading to a set of polyominoes. The general shape 
of the polyominoes is illustrated in Fig. [TU] It consists of a spike sent from the 
original square to the neighbor at distance n. 

As shown on the picture, communication bumps can be put on the spike, 
while dents are located inside the initial square. For matching the glues between 
one tile and its neighbor (i, 0), the bump will cross i—1 other spikes, we will see 
later how many layers of microtiles it represents. Indeed, the essential part of 
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Figure 10: General shape simulating a tile in arbitrary linear neighborhood. In this example, 
the neighborhood is N 3 = {(-3,0), (-2,0), (-1,0), (1,0), (2,0), (3,0)}. The shape is drawn 
in thin black, and the vertical thick lines are the communication places with the neighbor 
whose abscissa is the number indicated above it. These neighbors are drawn in gray. 
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Figure 11: Detail of the polyomino used for arbitrary linear-neighborhood zippers. Here, 
n = 2, I = 7, s = 28 and h = 6. The dark gray contour is filled by microtiles, the light gray 
area is filled for clarity and does not necessarily contain microtiles. 



the simulation is the discretization of this shape into a polyomino of microtiles. 
A final polyomino is represented in Fig. [TT1 

We define the following variables. The height of the polyomino is s microtiles, 
it is also the width of the initial tile and therefore the scale of the final poly- 
ribbon. The vertical space between the top of the polyomino and the beginning 
of the spike is denoted h, and since the spike is centered, the space between the 
end of the spike and the bottom is also h. Finally, the spike is a succession of 
horizontal segments of I microtiles, thus the slope of the spike is ±1/1. Then, 
the following constraints have to be respected when constructing the polyomino. 

• The integer I has to be as large as necessary, leaving enough horizontal 
space for crossings, as suggested by Lemma 13.31 Similarly, h needs to 
provide enough vertical space for the dents. Formally, I > a and h > (3, 
where a and /3 will be given later on. 

• The polyomino must be at least 3 layers wide everywhere, two for the 
inner and outer layers of the spike, and one for a potential junction of the 
path from the input microtile to the output microtile (as in Fig. | ^c)[ ). 
As a consequence, h > 3 and s > 3/ + 3, since the inner layer of the 
spike must go down by 3 before reaching s — 3. The first constraint can 
be removed (assuming j3 > 3), while the second one can be replaced by 
s > Al (provided I > 3, which is the case if a > 3), so that s can be exactly 
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divided in 4 horizontal segments everywhere on the spike. 
• The initial tile being a square, the height is s and can also be written 

h + 2[ns/l\ + h (the spike goes ns microtiles to the right at slope 1//). 

Therefore, after replacing both s by 42, we have h = 21 — An. 
Since h > j3, because of the last equation I has to be greater than 2n + (3/2. For- 
given n,a., ft E N, a solution of the system is then 

I = |max(a,2n + /3/2)] 
s = 4/ 

h = 2l-An . 

It is quite obvious that translated copies of this polyomino tile the plane with 
no overlaps, allowing to replace a grid of tiles by a grid of polyominoe^|. 

We now give some details on how the communications take place. For a 
more convenient description, we split the polyomino into n + 1 horizontal parts 
of width s = Al, we denote them from left to right by part (which corresponds 
to the initial square tile) to part n (end of the spike). Each of these parts is 
divided into 4 horizontal segments of length I, denoted segment 1 to segment 4. 
As suggested in Fig. [TU1 the bumps and bridges are located on the horizontal 
segments on the top of the spike in parts 1 to n, while the corresponding dents 
are on the horizontal steps on the top of the hole in part 0. For every < i < n, 
the glue ti t o of the initial tile is allocated some space in every part, on one of 
the four segments. A simple way to do this is to allocate glue t^o to segment 
(i — 1 mod 4), hence each segment is used for [n/A\ or |~n/4] glues. This space 
is then used in part for a dent, in parts 1 to i — 1 for bridges, and in part i 
for the bump. 

Moreover, there are i — 1 spikes to cross by the bump encoding tifi. Each 
spike is 4 layers thick, hence there are at most A(n — 1) layers to cross. This 
number is fixed, so we can apply Lemma 13.31 with k = 4(n — 1): each glue 
needs a horizontal space of (8(n — 1) + 4) (width of a dent) + (\X\ — l)(4(n — 
1) + 1) (spacing between different possible glues) microtiles. This gives a lower 
bound to I, which has to be greater than [n/4] x (\X\ (An — 3) + An — 1). After 
adding 6 microtiles to separate the dents from the sides of the polyomino, we 
define 

a = \n/A~\ x (\X\ {An - 3) + An - 1) + 6 , 

as the lower bound for I used previously. On the other hand, the depth of 
the dents is k + 1 = An — 3. Consequently, we have another constraint on h 
which must be greater than (3 = An (space for the biggest dent plus the three 
original layers). This means I greater than An, which is already the case because 



3 A formal proof of this fact could be made using the characterization of polyominoes tiling 
the plane given in [(J. Indeed, the contour word of our polyominoes (without bumps and 
dents, and filled to match the definition in Q) would prove them to be pseudo-squares. 
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I > a > An. Then, 

I = \n/A~\ x (\X\ (An - 3) + An - 1) + 6 
/i = 2Z - An . 

This ensures that our polyomino can be constructed, using an appropriate 
set of microtiles which will generate the X^-ribbon we described. The last step 
of the construction is the positioning of the input and output microtiles. We 
can choose to place them at top-left position (path coming from or going to the 
West), top-right (path from or to the East), middle of the top side (path from 
or to the North) , middle of the bottom side (path from or to the South) . Since 
s is even, the middle of a side is chosen after [s/2\ microtiles. Then the inner 
layer we preserved can be used for joining the input and the output microtiles 
easily; for example, starting from the input microtile, one can draw the contour 
of the path from the left, just before reaching the output microtile, the path 
goes one layer inside and goes back to the input microtile where it draws the 
contour from the right (see Fig. [T2] for examples). 



Figure 12: The three different paths when the input direction is West. 

Finally, putting the appropriate polyominoes one after the other generates a 
poly-ribbon which is in fact a Tjj-ribbon. As it was shown before, any X^-ribbon 
can be transformed into a poly- ribbon with \<jn horn which one can restore the 
original (T, iV„)-zipper using ip^ 1 . The poly- ribbon does not overlap if and only 
if the glues from the (T, 7V n )-zipper match everywhere, hence ip is a simulation. 

□ 

3.3. Complexity of the Linear- Neighborhood Construction 

We now give results on the "size" of this simulation, which underline the 
fact that the generated polyominoes can be very complex. All the results refer 
to the construction described in the proof of Theorem 13.41 

Lemma 3.5. A polyomino used to simulate a tile of a (T, N n )-zipper contains 
B = 0(n 2 ) bridges. 

Proof. When n = 2, there is one bridge between the neighbors at (—1,0) and 
at (1,0); when n = 3, there are in addition two bridges between neighbors 
(—2,0) and (1,0), and (—1,0) and (2,0). In general, there are n — 1 bridges 
between (i — n, 0) and (i,0) for 1 < i < n — 1, plus n — 2 bridges between 
(i — n — 1, 0) and (i + n — 1, 0) for 1 < i < n — 2, and so on. Hence there are 
B = Ym=i * = \ n i n ~ 1) bridges. □ 
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Lemma 3.6. A polyomino used to simulate a tile of a (T, N n )-zipper is consti- 
tuted by 0(\X\ n 3 ) microtiles. 

Proof. Drawing the contour of a shape requires: 

• s microtiles for the 2 horizontal segments in part 1; 

• Ah microtiles for all 4 vertical portions; 

• A(ns + An) microtiles for the 2 spikes (ns for the horizontal segments, An 
for the steps down and up); 

• at most s + 2h + 2(ns + An) + s/2 microtiles for joining input and output 
microtiles (worst case when joining left to bottom) ; 

• 2n bumps and dents of height at most 0(n) microtiles; 

• B = 0(n 2 ) (Lemma 13. 5[) bridges of height at most 0(n) microtiles, each 
one at most 3 layers thick. 

Since s — Al and we can choose I = \a] = 0(\X\ n 2 ), after summing all of the 
above we obtain the result. □ 

Lemma 3.7. Every tile of the initial (T, N n )- zipper is simulated by 0(\X\n 2 ) 
different polyominoes. 

Proof. For each tile t £ T, there are 4 (number of input positions) x 3 (number of 
output positions) x B \X\ (number of different possible bridges) different paths, 
where B is the number of bridges on the path. Since B = 0(n 2 ) (Lemma 13. 5p , 
there are C(|X| n 2 ) different polyominoes for one tile when n > 2. When n = 1, 
there are no bridges and there are only 12 different polyominoes. □ 

Proposition 3.8. In our construction, a T^-ribbon simulating a (T, N n )-zipper 
needs |T M | = 0(\T\ ■ \X\ 2 n 5 ) microtiles and = 0(\T\ ■ \X\ 2 n 5 ) glues. 

Proof. A (T, A r n )-zipper can use at most \T\ tiles, according to our construc- 
tion each of them is simulated by C(|X| n 2 ) different polyominoes (Lemma l3.7p 
constituted by C(|X| n 3 ) unique microtiles (Lemma l3.6p . Hence the simulation 
needs |T M | = 0(\T\ ■ \X\ 2 n 5 ) different microtiles. 

Since a polyomino is a ribbon which has to be uniquely built, microtiles 
(except for the input and output ones) have 2 glues which can be found only on 
one other microtile. Therefore each of these microtiles introduce a new glue, and 
reuse another: there are at least \X^\ = C(|T M |) = 0(\T\ ■ \X\ 2 n 5 ) glues. The 
other two glues are null\ and nulh, and the input and output glues microtiles 
use glues from X, which does not change the order of magnitude of □ 

3.4-. Simulating Arbitrary Neighborhoods 

In this section, we prove the first of the two main results of this paper, namely 
that zippers in any neighborhood can be simulated by ribbons (Corollarv l3.10|) . 
First, note that in a similar way to Remark 13.41 any neighborhood N can be 
replaced by an equivalent rectangular neighborhood N m ^ n = £ Z 2 | < 

\i\ < m, < \j\ < n and (i,j) ^ (0,0)} containing N. 
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Theorem 3.9. Let T be a tile system in rectangular neighborhood N m ^ n , with 
set of glues X . There exist a simulation ip and a tile system in neighborhood 
N v n, with set of glues X^, such that any (T,N m ^ n )- zipper can be simulated by 
a T^-ribbon. 

Proof. The key of the proof is the generalization of the ip simulation from the 
proof of Theorem 13.41 to rectangular neighborhoods. The idea is to have a 
vertical succession of n + 1 spikes of length m, each of them being a "sheath" 
for the next one (Fig. [13]) . 




Figure 13: Shape simulating a tile denned in a rectangular neighborhood N m>n of size (2m + 
1) X (2n + 1) (here m = 3 and n = 2). The initial tile is the light gray square, while the path 
is filled with darker gray. 

The bumps used for communications are put all along the spikes, at the 
places indicated on the picture. The difficulty is to find places for the dents, so 
that they can be contained in a place where enough space can be reserved. This 
is achieved as follows. 

• North (from (0, j), for < j < n) and West (from (i, 0), for — m < i < 0) 
communications are received in the area marked B on the picture. This 
works very similarly to the linear-neighborhood case. 

• Northwest communications (from (i,j), for — m < i < and < j < n) 
are received in A. Again, it is not difficult to see how the information 
crosses the layers. 

• Southwest communications (from (i,j), for — m < i < and — n < j < 0) 
are received in C. This is slightly more complicated to understand, since 
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these dents are not located inside the initial square. Putting the dent at 
the bottom of the polyomino allows to simulate the Southwest communi- 
cations by Northwest communications, which is then easily achieved by 
positioning bumps on the spike. 
The key point is that the areas marked A, B, and C (in dark gray on Fig. IT5]) 
are scalable, both vertically and horizontally, so they can be made as big as 
needed for the dents. Indeed, we can denote as previously by s the width and 
height of the initial tile, hence the scale of the poly-ribbon. Let x be the width 
of A, B, C and of all the other horizontal subdivisions of s. Since there are 
two of these blocks for each of the first n vertical spikes, plus one for the last 
spike, it holds that s — (2n + l)x. As in Theorem l3.41 we need the spike to be 3 
layers wide, hence the slope l/l of the spikes is defined by I such that x > 3/ + 3. 
Again it is possible to decide that x = 41, i.e., a block is made of four descending 
horizontal segments. Finally, let h be the left height of the A-B-C areas. The 
fact that the initial tile is a square is expressed by s = 2h + 2\ms/l\ . We have 
the following equations: 



Besides, to ensure enough space for the bumps and dents, we want to make sure 
that x and h are as big as necessary, i.e., x > a and h > ft (the exact values of a 
and f3 will be given later). Once solved, the system gives h = (2n + l)(2l — 4m). 
Since h should be greater than /3 and x greater than a, I needs to be greater than 
max(a/4, (3 /(An + 2) + 2m). Thus, a solution to the system which guarantees 
that the polyomino can be constructed is 



A last technical difficulty is the description of how the spikes shrink to fit 
into the previous one. The general way to do this is illustrated on Fig. 1141 which 
zooms on the rightmost part of the first spike of a polyomino. 

The outer spike does not shrink because it does not need to, since the shrink- 
ing will happen at the bottom of the initial tile (see Fig. [T3")) . This is not neces- 
sary, but it allows a tiling of the plane without any holes between polyominocs. 
The other spikes shrink by x microtiles on the left and on the right by going 
down 4 tiles (remember that the slope of the spikes is l/l — 4/x). Remark 
the slight asymmetry of the shrinking, which has to take place after reaching 
height s/2 on the left, and just before on the right. With all these conditions 
respected, it should be clear that vertically and horizontally translated copies 
of this polyomino tile the plane with no overlaps. 

It remains to determine the bounds a and /?. An application of Lemma 13.31 
to all the spikes and dents gives us the maximal size of bumps and dents. It is 



x = 41 

s = (2n + l)x 

s = 2h + 2\ms/l\ . 
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Figure 14: Illustration of how spikes shrink from width ax to (a — 2)x. The darkest gray 
corresponds to the spike of the initial tile, lighter grays represent the spike of farther North 
neighbors. 

obtained for the diagonal communication between a tile and its neighbor located 
at (to, —n), which crosses n+(m— l)(2n+l) spikes of thickness 4 layers, hence a 
total number of k — 4(2mn+m — n — l) layers. It follows that j3 = fc + 3 to ensure 
a free layer between the top of a dent and the upper side of the polyomino. The 
bound a is more complicated to express. As for Theorem 13.41 a communication 
bump-dent needs 2fc + 4 horizontal space, plus an extra (|X| — l)(fe+l) microtilcs 
for the different glues. The size x has to allow m + n dents in B, mn dents in A 
and C. Hence, after adding 6 microtiles for preserving the borders of the block, 
a = max(mn, m+n) x (2fc+4+(|X| — l)(fc+l))+6, with k — &(2mn+m—n—l). 

The rest of the proof (positioning and joining the input and output mi- 
crotiles, bijection between a (T, A min )-zipper and a set of unique poly-ribbons) 
is unchanged from the proof of Theorem 13.41 □ 

Corollary 3.10. Let T be a tile system in arbitrary neighborhood N , with set 
of glues X . There exist a simulation ip and a tile system in neighborhood 
N v n, with set of glues X^, such that any (T, N)-zipper can be simulated by a 
T^-ribbon. 

Remark 13. II can be extended to the general case, hence any Ijj-ribbon repre- 
sents a unique (T, A r )-zipper. Also note that a result similar to Proposition 13.81 
could be stated, but it would be more complex and of little interest since the 
number of tiles would be a lot bigger. 

To illustrate the achievability (and the complexity) of this construction, 
Fig. [15] gives a complete example of a polyomino for the simulation of a linear- 
neighborhood (T, A^-zipper, with \X\ = 2. 

The general case is by far more complicated and a detailed picture would be 
difficult to understand. We put in Fig. [16] the shape from Fig. [131 surrounded 
by 8 other shapes (colored alternatively in light and medium gray). The com- 
munication places between the central shape simulating the tile located at (0, 0) 
and the 8 other shapes are shown, one of them is detailed at the bottom of the 
picture. 
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Figure 15: Rotated polyomino simulating a tile of a (T, 7V2)-zippcr, with \X\ = 2, I = 23, 
L = 92, h = 38, input direction is West and output direction is South. The microtiles are 
located in the gray layer, the input and output directions arc indicated by the black arrows. 
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Figure 16: In dark gray is represented the shape from Fig. 1131 simulating tile r(0,0). It is 
surrounded by 8 other shapes corresponding to 8 neighboring tiles. Below is a detailed view of 
the circled communication place, with a bump crossing 4 bridges of width 4. Since the tile at 
position (1,2) is absent, one bridge is missing. The other bridges contain either 2 or 3 layers 
of microtiles, depending on how the input and output microtilcs arc joined. 
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4. Extension to Arbitrary-Neighborhood Tilings 

We now explain how the above construction can be modified, in order to 
simulate arbitrary-neighborhood tilings by von Neumann-neighborhood poly- 
tilings, to obtain our second main result fCorollary |4.2l) . Before that, we describe 
a simpler construction used in Q , which can be adapted to do such a simulation, 
but only in the case of total tilings. We finally study some properties of the 
tilings which are preserved by our simulation technique. 

4-1- Polyomino Construction 

The construction described in Sect. [3] transforms a tile into a polyomino of 
microtiles, keeping the path unchanged at the macro-level. Hence, a similar 
construction can be used to prove another result, which states that arbitrary- 
neighborhood tilings can be simulated by von Neumann-neighborhood tilings. 
As usual, we first prove the result for rectangular neighborhoods, and deduce 
the general case from that. 

Theorem 4.1. Let T be a tile system in rectangular neighborhood N m>n , with 
set of glues X . There exist a simulation tp and a tile system 7), in neighborhood 
N v n, with set of glues X^, such that any (T,N m>n ) -tiling can be simulated by a 
(T^A^jv) -tiling. 

Proof. Only a few modifications have to be done to the general polyomino repre- 
sented in Fig. [TBI in the construction of Theorem l3.9l to transform the T^-ribbon 
into a (T^, A t ,]v)-tiling. First, the inner layer used to join the input and output 
microtiles can be removed, allowing thinner spikes. 

The main difference is that now, all 4 glues of the microtiles have to match 
with their neighbors, if present. The glues towards the inside of the polyomino 
(i.e., the glues used for adjacent microtiles from the same polyomino) can be 
chosen as previously, uniquely for each tile of each polyomino. For the glues used 
to match with other polyominoes, things are slightly more complex, since they 
will be used to enforce the correct position of the adjacent polyomino, to avoid 
misalignment. These glues can be chosen as pairs (x, y) g N 2 . For North and 
East glues, (a;, y) represents the position of the microtile inside the polyomino, 
counting from the bottom-left part of the polyomino for example (hence, the 
microtile at the top- left of the polyomino will have glue (l,s(n + 1)) on the 
North side). For South and West glues, this pair should be the position of the 
microtile on the neighboring polyomino, to ensure proper match. Therefore, 
the West glue of the top- left microtile will be (s,s(n + 1)), since it has to 
match with the top-right tile of the polyomino adjacent on the right. Note that 
1 < x < s(m + 1) and 1 < y < s(n + 1) (see Fig. [13]), therefore the number 
of these glues is bounded, once T and N m ^ n are fixed. These glue do not need 
to encode any information from the original set of glues X, since this matching 
will be enforced by the shape of the polyominoes. Their role is only to align 
properly the polyominoes. 
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Then, each (T, iV TO) „)-tiling is associated with a set of (T M , iV,, ^-tilings by 
this construction. Conversely, any (T^, iV^j\r)-tiling r' can be turned into a poly- 
tiling t" using the tile system T M in neighborhood N v ?f, by discarding incomplete 
and misaligned polyominoes using a function \ v . Then, starting from the poly- 
tiling t", one can recover the initial tiles at their respective positions, with 
the help of the function (p^ 1 applied to the unique sets containing each of the 
complete polyominoes. By construction, this (T, iV)-tiling r is the only valid 
tiling such that Xv( T ') t P( T )i therefore ip is a simulation. □ 

Remark that instead of increasing the number of glues to force a correct 
alignment of the polyominoes, one could have worked on their shape, adding 
more bumps and dents as in [10( . We made the choice of a shape as simple as 
possible, even if it increases the number of glues used. 

As explained at the beginning of Sect. 13.41 any neighborhood can be replaced 
by an equivalent rectangular neighborhood, hence the following corollary. 

Corollary 4.2. Let T be a tile system in arbitrary neighborhood N, with set 
of glues X . There exist a simulation ip and a tile system in neighborhood 
N v n, with set of glues X^, such that any (T, N) -tiling can be simulated by a 
(Tfj,,N vN ) -tiling. 

4-2. Properties Preserved by the Simulation 

The main interest of this construction is that it guarantees that some prop- 
erties of the initial tiling are preserved, and still verified in the final poly-tiling. 
For example, unlike what happens with the simple construction explained in 
Sect. 12.31 a partial (T, AQ-tiling is simulated by a partial (T„, iV jy)- tiling. The 
converse is also true, as well as some other important properties like periodicity 
and convexity, as explained in this section. 

Partial tilings. Clearly, from our construction, a partial tiling is simulated by 
partial tilings. For the converse, remark that our construction creates "holes" 
in a poly-tiling, because the polyominoes we build are hollow. Then, provided 
we fill the polyominoes by new tiles with unique glues, we obtain the following 
immediate result. 

Theorem 4.3. Let T be a tile system in arbitrary neighborhood N , with set of 
glues X. There exists a simulation ip such that every (T, N) -tiling r is total if 
and only if all (T M , N v n) -tilings in tp( T ) are total. 

Note that when r is total, ip(r) is a singleton. Indeed, the polyominoes 
of ip(i) differ only by the position of the bridges. When the tiling is total, all 
bridges are constrained by the neighboring polyominoes, hence only one element 
of the set ip(t) is possible. The final poly-tiling is therefore uniquely defined. 

Remark J^.l. Similarly, remark that a tiling is finite [resp., connected] if and 
only if all the poly-tilings which simulate it are finite [resp., connected]. As a 
consequence, the simulation of polyominoes is achieved by polyominoes only. 
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Periodic tilings. Periodicity is an essential notion in the classical tiling theory [8j, 
[20| . it should be preserved by a meaningful simulation. The following result 
shows that this is the case with our construction. 

Theorem 4.4. Let T be a tile system in arbitrary neighborhood N, with set of 
glues X . There exists a simulation ip such that every (T, N) -tiling r is periodic 
if and only if at least one of the (T^, N v n) -tilings in ip(r) is periodic. 

Proof. Clearly, if a (T M , AT„jv)-tiling as constructed in Theorem 14. II is periodic, 
restricting it to a poly-tiling preserves its periodicity, and then the (T, iV)-tiling 
it simulates is also periodic. 

Conversely, if a (T, AT)-tiling is periodic, then one of the poly-tilings which 
simulate it according to Theorem 14. II has to be periodic. In the case of a total 
(T, iV)-tiling, this fact is obvious. In the case of a partial tiling, it suffices to 
choose the poly-tiling with the correct bridge constraints at the borders, so that 
the periodic pattern of the tiling repeats all over. □ 

Note that if the periods of a tiling r are p v ,Ph £ N+, the periods of the 
periodic poly-tiling in ip(r) are sp„, sph £ N+, where s is the scale of the poly- 
tiling. 

Convex tilings. The different convexity notions (line, column, or both) are im- 
portant when dealing with polyominoes [H, 0], which are a particular case of 
tilings. Our construction does not preserve the convexity of a tiling, but it is 
possible to modify it in order to preserve line convexity in the case of linear 
neighborhoods. 

Proposition 4.5. Let T be a tile system in arbitrary linear neighborhood N , 
with set of glues X . There exists a simulation ip such that every (T, N) -tiling r 
is line convex if and only if all (T^, N v n) -tilings in iJj(t) are line convex. 

Sketch of proof . First, polyominoes should be filled as in the proof of Theo- 
rem 14.31 Then, the construction needs to be modified again, for one reason 
visible on Fig. [FJ the bumps, dents and bridges create horizontal holes. 

The idea is to replace these "holes" by two sets of stairs, using again the 
position of the stairs to ensure the glue matching, as represented in Fig. [17] In 
this figure, the integer g is the position of the dent, as induced by some glue. 
In order for the top and bottom polyominoes to match, the bump has to be 
located at position g + k + 1, where k is the number of layers crossed. In the 
new construction, the glue matching is enforced by the fact that in the bottom 
polyomino, if g is too big then there is an overlap on the left stairs, and if g is 
too small then there is an overlap on the right stairs. 

The problem with this technical step is that crossing k layers implies that 
the spike goes down by 2(k + 1). If 7 is the number of communication places 
in a segment of length I, and S is the maximum number of layers which can be 
crossed, then each segment of length I in the spike would go down by 27(5+1) 
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Figure 17: Replacing bumps and dents by two stairs. Here, there are k = 2 layers to be 
crossed. 



instead of 1. This imposes new constraints on h and s, but the discretization 
system is still solvable, at the additional cost of increasing I by a factor 27(<5+l). 

Then, obviously, all polyominoes from the sets <p(t) are line convex. More- 
over, all polyominoes in any set tp(t) have the same height s and are assumed to 
have the same reference microtile, so that ip transforms a line of tiles into a block 
of microtiles of height s. This block is formed by the horizontal juxtaposition of 
polyominoes, which, by construction, do not leave holes between them (other- 
wise, there is a glue mismatch somewhere, leading to some overlap somewhere 
else, see Fig. HT]) . Thus, this block is made of s lines of microtileswithout holes, 
and is therefore line convex. Since a line convex tiling is a vertical superposition 
of lines, its image by ^ is a vertical superposition of line convex blocks, hence 
another line convex tiling. Finally, the converse is trivial, if a (T M , JV„jv)-tiling 
is line convex then it simulates a line convex (T, iV)-tiling. □ 

A similar result can be stated for column convex tilings, provided the initial 
tiling uses tiles in a vertical neighborhood N: if (i, j) 6 N then i = 0. The proof 
follows the same sketch as for Proposition 14.51 with everything rotated by 90°. 

Proposition 4.6. Let T be a tile system in arbitrary vertical neighborhood N , 
with set of glues X. There exists a simulation ip such that every (T, N) -tiling t 
is column convex if and only if all (T M , N v ^)-tilings in ^(r) are column convex. 

5. Conclusions and Perspectives 

In this paper, we proved that any zipper formed with tiles defined in an 
arbitrarily complex neighborhood can be simulated by a ribbon obtained by the 
catenation of simple polyominoes. Each of the new microtiles used for the simu- 
lation has, when placed on a ribbon, only two neighbors: its predecessor and its 
successor on the path. A similar construction can be achieved for tilings, in order 
to replace arbitrary neighborhoods by the simpler von Neumann neighborhood, 
and at the same time preserve some properties of the tiling. 

A few research directions can extend this work. For example, the simulation 
relies on the fact that zippers and ribbons are not self-crossing. Although this 
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is the standard way to define them, they could be generalized to other, self- 
crossing, notions. It would be interesting to see if similar results could be 
obtained in this case, using new constructions. Another interesting topic would 
be to study how this construction behaves in three-dimensions, since crossings 
might be avoided with spikes turning around each other. However, new problems 
occur and would have to be solved by original techniques. 

Another possible improvement is justified by the observation that our con- 
structions are not yet adapted to practical dynamical implementations, such as 
DNA self-assembly. Theoretically, our constructions could be used for practical 
purposes: for example, as mentioned in [21|, the construction of Theorem 13.11 
can be used to simulate a Turing machine at temperature 1 (i.e., DNA tiles at- 
tach to the growing assembly whenever one glue matches). But this produces a 
lot of "garbage" consisting of many blocked polyominoes, for a limited number 
of correct assemblies. Indeed, there are many dynamical scenarios wherein, for 
example, the self-assembly of a "bad" polyomino is started, that blocks any fur- 
ther aggregation. This issue could be solved by adding the possibility to recover 
from a wrong start, for example by allowing tiles to "unstick", as suggested 
in }28l [l|. Further modifications would then be necessary to guarantee that, 
also in this setting, our constructions can self-assemble fully and correctly in 
finite time. 

Also remark that in our constructions, we emphasized the regularity of the 
polyominoes, to provide easier general constructions and aggregation. This led 
however to a blow-up in the number of new microtiles necessary to simulate a 
zipper, potentially leading to more experimental issues. The next step would 
be to optimize the construction according to some of the following criteria: 
number of polyominoes associated with a tile, number of microtiles appearing 
in a polyomino, number of glues used to build the polyominoes, etc. 

An alternative solution to these two problems would be the use of staged 
self-assembly [l2| for experiments. This formalism would allow to add an initial 
stage to construct the polyominoes in separate bins, before mixing them in the 
final solution, preventing the formation of "bad" assemblies. Besides, staged 
self-assembly would dramatically decrease the number of required glues (hence 
the number of microtiles), by adding even more initial stages to produce the 
polyominoes. This would offer more control on the order in which microtiles 
attach, thus removing the necessity for unique glues. However, this would be 
achieved at the cost of increased control by an external operator during the self- 
assembly process. This trade-off between external control and tile complexity 
is often encountered in DNA self-assembly (see, e.g., HI 12 1), and is currently 
being investigated by the authors. 

The authors would like to thank Shinnosuke Seki for his useful suggestion 
simplifying the construction in Fig. 1131 and Leonard Adleman, Jarkko Kari, and 
Erik Winfree for discussions. 
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